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Abstract 

On a smooth bounded domain 17 C we consider the Schrodinger operators — A — V^, with V being 
either the critical borderUne potential V{x) — {N — 2)^/4 |a;|~^ or V{x) = (1/4) dist(x, 917)^^, under 
Dirichlet boundary conditions. In this work we obtain sharp two-sided estimates on the corresponding 
heat kernels. To this end we transform the Schrodinger operators into suitable degenerate operators, 
for which we prove a new parabolic Harnack inequality up to the boundary. To derive the Harnack 
inequality we have established a series of new inequalities such as improved Hardy, logarithmic Hardy 
Sobolev, Hardy-Moser and weighted Poincare. As a byproduct of our technique we are able to answer 
positively to a conjecture of E. B. Davies. 
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1 Introduction and main results 

Harnack inequalities have been extremely useful in the study of solutions of elliptic and parabolic equations, 
starting from the pioneering works of De Giorgi |D(Tj . Nash and Moser |Molj . |Mo2j . They are used 
to prove Holder continuity of solutions, strong maximum principles, Liouville properties, as well as sharp 
two-sided heat kernel estimates. In particular, we should mention the influential works of Aronson 
and Li and Yau |LYj where heat kernel estimates were obtained via parabolic Harnack inequalities. 

In fact, in certain cases, the parabolic Harnack inequality is equivalent to sharp two-sided heat kernel 
estimates. This is the case when dealing with second order uniformly elliptic operators in divergence form 
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on R-^, or more generally with weighted Laplacians on complete Riemannian manifolds; see the works 
of Fabes and Stroock |FSj . Grigoryan |G1| . and Saloff-Coste |SClj . This equivalence has been also used 
in order to get sharp two-sided estimates for Schrodinger operators in R^. For instance, the case of a 
potential that is regular and decays like at infinity was studied by Davies and Simon where 

pointwise upper bounds for the heat kernel were derived. The picture was later completed by Grigoryan 
|G2j where sharp two sided estimates were provided by means of a parabolic Harnack inequality. A recent 
survey on heat kernels on weighted manifolds can also be found in |(T2j . 

As it was shown in the works of Fabes, Kenig and Serapioni |FKSj . and Chiarenza and Serapioni 
|('Sj . parabolic Harnack inequalities follow after establishing Poincare and Sobolev inequalities as well as 
a doubling volume growth condition. Moreover, on complete Riemannian manifolds parabolic Harnack 
inequalities are equivalent to Poincare inequality and a doubling volume growth condition as explained 



by Grigoryan and Saloff-Coste in [(jSC'j . |S(]2j . 

Since the work of Baras and Goldstein |BGj , the existence or nonexistence of solutions to the partial 
differential equation 

ut = Au + Vu, (1.1) 

with a potential V involving the inverse square of the distance function have been widely investigated. See 
|BGj . Brezis and Vazquez |BVj . Cabre and Martel |CMj . as well as Vazquez and Zuazua |VZj . for the case 
V{x) = c|x|~^ and |CMj for the case V{x) = cd~'^{x) on a bounded domain 17, where d{x) = dist(x, dO,). 

Concerning the case where V{x) = c|x|~^ with c < (A^ — 2)^/4 , sharp two-sided heat kernel estimates 
have been obtained in R'^, see |MTlj . |MT2j where the approach of (GSCj on complete Riemannian 
manifolds has been used, after a suitable transformation; see also |MSj for a different method. 

On the other hand few results are known in the case of incomplete Riemannian manifolds, as it is for 
example the case of bounded domains in R^. To our knowledge the only sharp two sided estimates in 
this case, concern the standard Dirichlet Laplacian on a smooth bounded domain C R^, first studied 
by Davies and Simon in |Dlj , |D2j , |DSlj , and recently completed by Zhang [Zj . We note that in the case 
of a bounded domain, the asymptotic of the heat kernel is different for small time than it is for large time. 
In fact, for the heat kernel of the standard Dirichlet Laplacian and for two positive constants 

Ci ^ C2, we have for small time 

^ . d(x)d(y)} N ix-i/p N ^ . f. dix)diy)\ n ^, \^-y\^ . . 

Cimm<{ 1, ^ \ U'-e'^'—r- < hD{t,x,y) < C2 mm <^ 1, ^ U'-e"^^—^, (1.2) 



whereas for large time 



Ci dix) d{y) e-^i* < hDit,x,y) < C2 d{x) d{y) e'^^*, (1.3) 



for all x,y £ i}; here Ai is the first Dirichlet eigenvalue. 

In this work, our main interest is in obtaining sharp two-sided estimates for the heat kernel of the 
Schrodinger operator —A — V under Dirichlet boundary conditions, on a smooth bounded domain O C R'^ 
for the following critical borderline potentials: V{x) = {{N — 2)^/4)|x|~^ or V{x) = {l/4)d~'^{x). 

Throughout this work $7 is a bounded domain of R^ containing the origin and d{x) = dist(x, dO,). 
We first consider, for A > 3, the case V{x) = ((A^ — 2)^/4)|x|^^, x G and we formally define the operator 
K by Ku = —Au — ^^^J^-} u , u\qq = . More precisely, the Schrodinger operator K is defined in L^(J7) 
as the generator of the symmetric form 



namely, if 



f f (A -2)2 \ 

}C[ui,U2] := / VniVn2 -j— 75 — U1U2 dx , 

Jn V J 

D{K) := j-u G H{Q) : -An - \ G L2(0)| 
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(N - 2)2 

Ku := -Au - ^ — no^u for any u £ D(K) , (1.4) 
4 xr 



where H{Q) denotes the closure of C^(r2) in the norm 



u\\Hin) ■.= \ ii |Vn|2 - '-^r^u' ]dx\ . (1.5) 



(7Vj-2^2 

4|a;|2 

Let us recah that H{^1) C Wq''^{Q) for any I < q < 2, due to the results in Subsection 4.1 of |VZj . 

It follows, using Hardy inequality, that is a nonnegative self-adjoint operator on L'^{Q) such that 
for every t > 0, e~^* has an integral kernel, that is, e~^*tio(x) := f^k{t,x,y)uo{y)dy where k{t,x,y) is 
the heat kernel of K. The first Dirichlet eigenvalue of K can be defined by 

In (IV^I' - ^'^) ^ 
Ai := mf ^ J. — 5- , (1.6) 

with Ai > 0, due to [BV]. Moreover there exists a positive function cpi G H{Q) satisfying 

{N - 2)2 

— Awi — — ifi = Ai^i, in ipi = 0, on dQ, 

see for example Davila and Dupaigne |DDj . 

We then have the following sharp two-sided heat kernel estimate on K for small time 

Theorem 1.1 Let C R^, > 3, be a smooth bounded domain containing the origin. Then there exist 
positive constants Ci,C2, with Ci < C2, and T > depending on Q, such that 

^ . f/i I ^- -I , 1-^ N-2 d(x)d(y)), , 2-N N ^ \^-y?- 
Cimm|(|a;| + ^ft)~{\y\ + Vt)~, ^ ' ' Y^y\~ t~~ e~^^^^ < 

^ f/i , ^. N-2 ,. , ^, N-2 d(x)d(y)"\, . 2-N N \x~y\'^ 

< k{t,x,y) < C2mm|(|x| + Vt) — {\y\ + Vt) — , --^^-^ j\xy\~t~- e~'^^~^ , 
for all x,y £ il. and < t < T . 
Concerning the large time asymptotic we have: 

Theorem 1.2 Let C R^, > 3, be a smooth bounded domain containing the origin. Then there exist 
two positive constants Ci,C2, with Ci < C2, such that 

Ci d{x) d{y) |xy|^e~^i* < k{t,x,y) < C2 d{x) d{y) \xy\^ e'^^\ 

for all x,y £ il. and t > large enough; here Ai is defined in M.b}) . 

To prove the above Theorem ll.2l we have shown a new improved Hardy inequality which is of independent 
interest; see Theorem 13.21 

We next consider the case where the Schrodinger operator H has a potential with critical borderline 
singularity at the boundary Hu = —An — j^^j^u , u\gQ = 0; here N > 2 and 17 is a convex domain. 
More precisely, the Schrodinger operator H is defined in L'^(Q) as the generator of the symmetric form 

U1U2 ] dx , 



n[u„U2] := j^{vu,Vu2-j^^u,U2^ 
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namely if 



DiH) := |n E : -Au - G L^Q)^ , 

Hu := —Au — u for any u G D{H) , (1-7) 

where denotes the closure of Cq°{^}) in the norm 

u ^ WuW^^n) := { I (|Vn|^ - dx^ ' 

Let us recah that W{n) C W^'''{n) for any 1 < g < 2, due to Theorem B in |BFTT]. 

Then, due to Hardy inequality, H is a nonnegative self-adjoint operator on L^(r2) such that for every 
t > 0, e~^* has an integral kernel, that is, e~^*uo(x) := f^h{t, x,y)uo{y)dy; here h{t,x,y) denotes the 
heat kernel of H. The first Dirichlet eigenvalue of H is defined by 



Ai := mf — 5- . (li 

o^<pec^{n) Jn^ dx 



It is known that Ai > —00 for any bounded domain Q,, and Ai > if is convex, see |BMj . Moreover 
there exists a positive function ipi £ W{^1) satisfying 

1 

-Aifi - j-j^j-^ipi = Xiifi, m 9^1 = 0, on dQ; 

see for example |DDj . 

We then have the following sharp two-sided heat kernel estimate on H for small time 

Theorem 1.3 Let Q C R^, N > 2, be a smooth bounded and convex domain. Then there exist positive 
constants Ci, C2, with Ci < C2, and T > depending on such that 
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^ I d2(x)d2(y)\ N ^, \x-v\^ ^ ^ 1 d2(x)d2(y)\ N \x-y, 

Cimin<^ 1, — ^ \ U-^e-^^^- < h{t,x,y) < Camin j 1, — ^ \ \t'^^^^^^^ ^ 
for all X, y £ 0, and < t < T . 

We next complement this with the large time behavior: 

Theorem 1.4 Let Q C R^, N > 2, be a smooth bounded and convex domain. Then there exist two 
positive constants Ci,C2, with Ci < C2, such that 

Ci d^{x) d^{y) e^^i* < h{t,x,y) < C2 d^x) d^{y) e"^!*, 

for all X, y £ Q and t > large enough; here Ai is defined in U.^) . 

The two-sided estimates in Theorems 11.11 and 11.31 are obtained as a consequence of a new parabolic 
Harnack inequality up to the boundary, for a suitable degenerate elliptic operator. Let us present a model 
operator in this direction. For this we consider classical solutions of 

vt = ^^dzv{d'^{y)Vv), (1.9) 



(actually solutions are considered as weak solutions, for the precise formulation we refer to Definition 12.91 
with A = there, note that due to elliptic regularity, any solution is smooth away from the boundary of 

n). 

Then, the following Harnack inequality holds true: 
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Theorem 1.5 (Parabolic Harnack inequality up to the boundary). Let N > 2, a > 1 and 

Q, C be a smooth bounded domain. Then, there exist positive constants Ch and R = R{^) such that 
for xGil, 0<r<i? and for any positive solution v{y, t) of !ll.y\) in {B{x, r) nO,} x (0, r^), the following 
estimate holds true 



ess sup 



{y,i)e{s(.,i)nf7}x(4,^) ^(2/'*)^^^^ e«^^^/(y,t)e{B(x,i)nn}x(|.^.2) ^V^^) ■ (1-10) 

Here B{x, r) denotes roughly speaking an N dimensional cube centered at x and having size r (see Defini- 
tion 12.1(1 . The estimate of Theorem 11.31 is a consequence of Theorem 11.51 and corresponds to the extreme 
value a = 1. We refer to Theorem 12 . 1 ( )l for a more general result that leads to Theorem ll.il 

The existence of a uniform upper bound on the size of the admissible "balls" denoted by R{Q), 
is necessary, because otherwise the nonexistence of an upper bound would imply two-sided heat kernel 
estimates that are the same for small time and large time, which is not the case at least for a = 1, due to 
Theorems 11.31 and 11.41 

The restriction on a in Theorem 11.51 is sharp, since in the weakly degenerate case, where < a < 1, 
even the elliptic Harnack fails. Indeed, let := B{0, 1), then v{y) := is a positive solution 

of div{d^{y)^v) = for 1/2 < \y\ < 1, with v{l) = 0. The natural analogue of Theorem ll.Sl in the weakly 
degenerate case, that is < a < 1, is a Harnack inequality for the ratio of any two positive solutions; in 
the elliptic case this is done by a probabilistic approach in 

To derive heat kernel estimates we define the operator L := — ^J^^^ div{d" (x) V) in L'^{Vl,d°^{x) dx) as 
the generator of the symmetric form 



C[vi,V2\ := / (i°(x)Vf iVu2 dx 

namely 



D{L) := |u G Hlin^d'^ix) dx) : -^^^div{d'^{x)Vv) G L2(0,d°(x) dx)^ , 

Lv := — --^div(d°'(x)Vv) for any v G D(L) , (1.11) 
a"(x) 

where //g (fi, (i"(x) dx) denotes the closure of C^(J7) in the norm 

y^\\y\\^,:=n d"(x) (|V7;|2 + ^2) dxV . (1.12) 



We should emphasize that for a > 1, one has H^{n,d'^{x) dx) = H^{n,d°'{x) dx), see Theorem ITTTl 
Let us note that L is a nonnegative self-adjoint operator on ^^(rj, d"{y)dy) such that for every t > 0, 

e~^* has a integral kernel, that is e~^^vo{x) := f^l{t, x,y)vo{y)d'^{y)dy; the existence of the heat kernel 

l{t, X, y) can be proved arguing as in (DSlj . 

Then we obtain the following sharp two-sided estimate for the heat kernel generated by L. 

Theorem 1.6 Let a > 1, N > 2 and C be a smooth bounded domain. Then there exist positive 
constants Ci,C2, with Ci < C2, and T > depending on Q such that 

^ . f 1 1 ] -E. 1"^-^!^ ,/ \ ^ ■ f 1 1 1 _iv \^-y\^ 

C;imm<^ — , W 26^2 t < l[t,x,y) < C;2mm<^ — , W 2 e t , 

Lt2 d2[x)d2[y)) \^t2 d2[x)d2[y)) 

for all x,y £ il. and < t < T . 
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We also suppose that 



So far we have considered special potentials V. However as we shall see next we can obtain much 
more general results. For instance we consider the operator E := —A — V where the potential V is such 
that 

V{x)=Vi{x) + V2{x) , (1.13) 

where 

l^i(a^)l<^^, V2{x)eLPin), p>^ . (1.14) 

L (|V^|2 - Vif'^) dx , ^ 

Ai:= inf ' >0, 1.15 

and that to Ai there corresponds a positive eigenfunction ipi satisfying for all x G 1^ the following estimate, 

cid"2 (x) < ipi{x) < C2(i2"(x) , for some a > 1 (1-16) 

and for ci, C2 two positive constants. 

Then as before it can be shown that £^ is a well defined nonnegative self-adjoint operator on L'^{Q) 
such that for every t > 0, e~^^ has a integral kernel, that is e~^^uo{x) := e(t, x, y)uo(y)dy. We consider 
positive solutions of 

ut = -Eu ; (1.17) 

then our first result reads 

Corollary 1.7 For N > 2, let Q C be a smooth bounded domain. Suppose that lll.l!^) . \l.i4\ l, il.l5]) 

and il.16]) are satisfied. Then, there exist positive constants Ch o-nd R = R{Vt) such that for x ^ Vl, 
< r < R and for any positive solution u{y,t) of jl-lTj ) in {B{x,r) PiO,} x (0,r^) we have the estimate 

'''P(y,t)e{B{.,^,)nn}xi^4) u{y,t)d''l (y) < Ch ess ^n/(^_,)g|^(^,.)n^|^( 3,2^,2) u{y,t)d-^{y) . 

Our result in the case a = 2 is basically the local comparison principle by Fabes, Garofalo and Salsa 
|FGSj in the case of Schrodinger operator (see Remark 12.161 that covers the uniformly elliptic case). 

As usual from the parabolic Harnack inequality the following sharp two-sided estimate for the heat 
kernel e{t,x,y), corresponding to the operator E, can be easily deduced: 

Corollary 1.8 For N > 2, letQc be a smooth bounded domain. Suppose that m.l!J\) . \1.14\ ), HI .15]) 

and il.16]) are satisfied. Then there exist positive constants Ci,C2, with Ci < C2, and T > depending 
on such that for any x, y G 17 and < t <T 



^ d2(x)d2(y) N \^~y\^ , s ^ . L d^(x)d^(y) n \^-y\^ 
1, — ^ L \t--e~^^—r- < e{t,x,y) < C2 mm <^ 1, — ^ W^Te"^^^^, 

whereas for t > T we have 

Ci dt(x)dt(y)e-^i* < e{t,x,y) < C2 d^ {x)d^ {y)e-^^^ . 

As a byproduct of our method we can answer a conjecture by E. B. Davies (Conjecture 7 in |D2j ) in the 
case of the Schrodinger operator (see Section 4.4 for a more general case). For this let us introduce the 
Green function associated to E, that is 



GE{x,y) = / eit,x,y)dt , (1.18) 

then we have 
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Corollary 1.9 For N > 3, letQc be a smooth bounded domain. Suppose that M.ltl\) . \l.lJi\j , 111.15}) 

and il.ld]) are satisfied. Then there exist two positive constants Ci,C2, with Ci < C2, such that for any 
x,y £ Q, 

Davies conjecture corresponds to our result in the case a = 2, we should note however that other 
values of a > 1 are possible. 

The structure of the paper is as follows. In Section 2 we prove the new parabolic Harnack inequality 
up to the boundary for a doubly degenerate elliptic operator as well as the two sided heat kernel estimates 
that can be deduced from it. In Section 3 we present the proof of the above mentioned results concerning 
the Schrodinger potential having critical singularity at the origin, while Section 4 treats the case of the 
Schrodinger operator having critical singularity on the boundary. 

Acknowledgment This work was largely done whilst the second author was visiting the University 
of Crete and FORTH in Heraklion, the hospitality of which is acknowledged. This research has been 
partially supported by the RTN European network Fronts-Singularities, HPRN-CT-2002-00274. 
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2 Parabolic Harnack inequality up to the boundary for degenerate 
operators 

In this section we prove a new parabolic Harnack inequality up to the boundary for the doubly degenerate 
elliptic operator in divergence form 

L'a ■■= -p^d..(|x|^d"(x)V), (2.1) 

for any a > 1 and AG [2 — A^, 0]. To this end we will use the Moser iteration technique. According to the 
approach presented in |GSCj as well as in |CSj the key estimates one needs are the doubling volume-growth 
condition (see Corollarv 12. 4|) . the local weighted Poincare inequality (see Theorem 12. 5() as well as a local 
weighted Moser inequality. In fact we will establish two local weighted Moser inequalities, one will be 
used near the boundary (Theorem 12. Hj) and the other one away from the boundary (Theorem 12.13)1 . 

Then, similarly as in |GSCj . we deduce from it a sharp two-sided heat kernel estimate. To this end a 
sharp volume estimate is needed (see Lemma l2. 2)) . 

In the sequel we will use the following local representation of the boundary of fi. There exist a finite 
number m of coordinate systems {y[.,yiN), ui = {yn, ■ ■ ■ ,yiN-i) and the same number of functions Oj = 
ai{y[) defined on the closures of the {N — 1) dimensional cubes Aj := {y'^ : \yij\ < 5 for j = 1, • • • , — 1}, 
i £ {I,-"" ;"^} so that for each point x £ there is at least one i such that x = {x[,ai{x'j)). The 
functions Oj satisfy the Lipschitz condition on Aj with a constant ^ > that is 

Wiiy'i) - ai{zi)\ <A\y[-z[\ 

for y\^z\ € Aj; moreover there exists a positive number /? < 1 such that the set Bi defined for any 
i G {1, • • • ,m} by the relation Bi = {{y\,yit^) : y[ £ Ai,ai{yl)-(3 < yiN < ai{y'^) + l3} satisfy Ui = BiPiO, = 
{{y'i,yiN) ■■ y'i e Aj,aj(y-) - 13 < ytN < ai{y[)] and Fj = Sj n 9^7 = {{y[,yiN) ■ y'i e Aj,yjAr = aj(y^)}. 
Finally let us observe that for any y £ Ui one can make use of the following estimate on the distance 
function (1 -|- A)~^[ai{y[) — yny) < d{y) < {ai{y[) — yiN) (see Corollary 4.8 in ^ for details) 

Let us fix from now on a constant 7 £ (1, 2) and let us define the "balls" we will use in Moser iteration 
technique. Roughly speaking they will be Euclidean balls if they stay away from the boundary and they 
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will be dimensional "deformed cubes" , following the geometry of the boundary, if they are close enough 
to the boundary or even if they intersect it. More precisely we have 

Definition 2.1 (i) For any x G and for any < r < l3 we define the "ball" centered at x and having 
radius r as follows. B{x,r) = B(x,r) the Euclidean ball centered at x and having radius r if d{x) > "fr, 
while B{x,r) = {{y'i,yiN) ■ Wi - x'^l < r, ai{y[) - r - d{x) < yuq < ai{y[) + r - d{x)} if d{x) < -/r, where 
i G {1, • • • ,m} is uniquely defined by the point x G d^l such that \x — x\ = d{x), that is by the projection 
of the center x onto dQ. (ii) We also define by V{x,r) := Jig(^^ r)nn \y\'^d"'{y)dy the volume of the "ball" 
centered at x and having radius r. 

We first derive a sharp volume estimate. 

Lemma 2.2 Let a > 0, N > 2, X ^ (— A^, 0] and C be a smooth bounded domain containing the 
origin. Then there exist positive constants ci,C2 and f3 such that for any x £ Q and < r < f3, we have 

cimax{d"(x)(|x| +r)^,r"}r^ < V{x,r) < C2 max{(i"(a;)(|x| +r)^,r°}r^ . 

To this end we make use of the following Lemma which can be proved as in |MT2j . 

Lemma 2.3 Let N > 2, X (— A^, 0] and J7 C be a smooth bounded domain containing the origin. 
Then there exist two positive constants di,d2 such that for any x £ 0,, we have 

di r^{\x\ + r)^ < [ \y\^dy < r^(|x| + r)^ . (2.2) 

JB{x,r) 

Let us accept (|2.2|) at the moment and let us prove the sharp volume estimate. 

Proof of Lemma \2.^ ' Let us first consider the case where d{x) > jr. Then B{x,r) = B{x,r) C il. 
Due to the fact that any y G B{x,r) satisfies 

d{x) < d{x) - r< d{y) < d{x) + r< ( ) d{x) (2.3) 



7 / V 7 



the claim easily follows making use of Lemma 12.31 with C2 > d2 ( -^-^ ) and ci < d^ ' !- 
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Let us now consider the case where d{x) < jr and let us denote by Li,L2 two positive constants 
such that B(0,Li) C C B(0,L2) (note that they exist by assumption on 17). Then any y G B{x,r) D 
satisfies the following estimate Li — (7 + l)/3 < \y\ < L2. Indeed, if on the contrary \y\ < Li — (7 + l)/3, 
then by definition of Li we would have d{y) > (7 + 1)/?, and this contradicts our assumption. In fact one 
obviously has d{y) < d{x) + r < (7 + l)r < (7 + 1)/3 and it is not restrictive to suppose from the beginning 
that the parameter /? in the local representation of the boundary of Q satisfies f3 < ^1(7 + 1)~^. As a 
consequence we have: 

yy£B{x,r)nn d3<\y\^<d4. (2.4) 

here ^3 := L^ and ^4 := (Li - (7 + l)l3)^. 
Then for some i G {1, • • • , m}, we have 

V{x,r)= \y\>^d'^{y)dy ^ / \y\^{aM) - y^NTdyiNdy[ . 

JB{x,r)nn J \y'i-x'i\<r Jai{y'.)-r-d{x) 

From now on we omit the subscript i for convenience. Indeed we have 

/> i-mm{a{y'),a{y')+r-d{x)} p 

V{x, r)< / \y\\a{y') - ynTdyNdy' < ^4(7 + l)r{d{x) + r)" / dy' < 

J \y' —x'\<.r J a{y')—T—d(x) J \y' —x'\<r 
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< d4(7 + l)"+V"+^+^-^u;^_i = ^4(7 + l)"+V"+^u;^_i . 

On the other hand 

Vix, r) > (1 + ^)-" / / \y\\a{y') - yNTdyNdy' > 

J \y' —x'\<r J a{y')~i — d{x) 



r ra.{y )-rcy~l) 

>d3(l+^)-"/ / (r(7 - I)rd2/W > 

J\y'-x'\<r Ja{y')-r-d{x) 

>ds{l + A)-"r^+\2 - 7)(7 - 1)" / dy' = d^{l + A)-"r^+'+^-\2 - - l)"^iv-i 

J \y' —x'\<r 



\y'- 

= ds{l + ^)-°r"+^(2 - 7)(7 - ■ 
Here a; at denotes the standard volume of the Euchdean unit ball in R^. Thus the result follows with ci := 
minjdi (^^y ,d3{l + A)-^{2 - -f)i-f - iroJN-i} and C2 := max {da (2±i)" , ^4(7 + 1)"+ Wi}- 

// 

Let us now prove estimate 1)2. 2() , which is taken from |MT2j , we give here the details for the convenience 
of the reader 

Proof of Lem.ma \2.c\ (i) Observe that 

\y\^dy = r^+^ / \w + z\^dz , 



B{x,r) JBiO,!) 



where w := -. Then 1)2.2^ reads 



di{\w\ + l)^< \w + z\^dz<d2{\w\ + l)^ . (2.5) 

Jb{o,i) 



Since \z\ < 1, there holds 

liul — 1 < \w\ — \z\ < \w + z\ < \w\ + \z\ < \w\ + 1 



whence 



uJn(\w\ + iY < [ \w + z\^dz < un(\w\ -'^Y (2-6) 
^ ^ Jb(oa) ^ ' 



?(o,i) 

Comparing (|2.5|) and (|2.6|) . we immediately obtain the lower bound in (|2.2|) with d\ := uj^. (ii) To prove 

M, (b) M 



the upper bound in (|2.2() . observe that three cases are possible: (a) r < (b) ■'y < r < 3|x| and (c) 



r > 3|x|. In case (a) the claim follows from the right-hand inequality in (|2.6() . if we exhibit 1^2 > such 
that 

ujn{\w\ -1)^ <d2{\w\ + l)^ 

for any \w\ > 2. It is easily seen that the function F{t) := (^f^r) (i > 1) is decreasing, thus F{t) < 

F{2) = (1)^ for any t >2. This proves the claim in this case for any d2 > uj^-^- 
To deal with cases (6) — (c), observe first that 

B{0,r) A + iV 

In case (6) we have 

B{x,r) <ZB{0,4\x\) , (2.7) 
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whence 



In case (c) there holds 



JB{x,r) A + iV A + iV 

i?(a;,r) csj'cy') , (2.9) 



thus 



Since 



B{x,r) 



/ y, \ X + N 



A + V 3 / 



we have 



^A+^ < ^"^(1^^ + ^)' for |x| < 2r 

in cases (b)-(c). Hence the conclusion follows with d2 ■= ^\ ■ 

^Prom Lemma 12.21 one can easily deduce the doubling property which reads as follows: 



// 



Corollary 2.4 Let a > 0, > 2, A G (— A^, 0] and C be a smooth bounded domain containing the 
origin. Then there exist positive constants Cd and j3 such that for any x G and < r < 13, we have 

V{x,2r) < CDV{x,r). 
Let us state now the local Poincare inequality. 

Theorem 2.5 (Local weighted Poincare inequality) Let a > 0, N > 2, X (— A^, 0] and Q C be 

a smooth bounded domain containing the origin. Then there exist positive constants Cp and [3 such that 
for any x and < r < (3, we have 

inf / \y\^d"{y)\f{y)-C\^dy <Cpr^ [ |2/|^d"(y)| V/I^dy , V / G CHBix,r)nn) . (2.11) 

^^^JB{x,r)nn JB{x,r)nn 

Proof: Let us first consider the case where d{x) > jr. Then B{x,r) = B{x,r) C 0. Due to 1)2. 3() the 
claim corresponds to Theorem 3.1 in |MT2j . We give here the details for the convenience of the reader. 

(i) As a consequence of the compact embedding of the space H^{B{0, 1), \y\''^dy) into L^(i?(0, 1), \y\^dy) 
(e.g. see [KOn we have that 



/ \f-fWdy<c[ \Vf\'\y\^dy, V / G C7i(S(0, 1)) ; 
Jb{o,i) -^5(0,1) 

here / := (^Jb(o,i) f{y)\y\^dy^ (/b(o,i) 1^1^^^) • Then by scaling p.ll|) follows when a; = 0. 

(ii) Let us now consider the case \x\ > 2r and let us define / := /^^q f{x + rz)dz. Then as in 



the proof of Lemma 12.31 we have (|x| + r)^ < |x + rz\^ < ^(|x| + r)^. Hence 



[ \f-f\^\y\^dy = r^ [ \f{x + rz)-f\^\x + rz\^dz<^{\x\+r)^[ \f(x + rz)-f\^dz< 
JB{x,r) JB{0,1) -3 -^5(0,1) 
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J.N f- N r If 

3 -'B(0,1) 3 Jb(0,1) 3 JB{x,r) 

Then ((TTT|) follows when >2r. 

(in) For a general x £ Q two cases are possible (a) < < |; (b) |x| > |. In case (a) there holds 

^^(O,0 Ci?(x,r) , 

thus from (i) we have 



inf / \fiy)-C\^\y\^dy < inf / \fiy)-C\^\y\^dy < ^r^ [ \Vf\^\y\^dy < Cpr^ [ |V/| 

This proves H2.11|) in case (a) since using a Whitney type covering and arguing as in |SC2j the integration 
set of the left hand side which from above is B (x, |) can be increased as to cover all B{x,r). 
In case (b) there holds \x\ > 2 (0; hence from (ii) 

inf / \f{y) - ^\^\y\^dy < ^r' [ \Vf\^\y\^dy < Cpr^ [ \Vf\^\y\^dy . 

«GR7ij(^.,|) 64 Jb{x,^) JB{x,r) 

This completes the proof in the case d{x) > ^r. 

Let us now consider the case where d{x) < 'jr. Then for some i £ {!,■ ■ ■ , m} we have 

d''iy)\fiy)-^\^dy< / / \fiy^,y^N)-^\''\y\\a,iy',)-y,Mrdymdy', 

B{x,r)nn J\y'i-x'i\<'r J ai{y[)-r~d{x) 

From now on we omit the subscript i for convenience. Let us perform then the following change of variables 
{y' 1 yN) iy', zn ■= o,{y') — Vn) and make use of H2.4|) : thus the above integral is less or equal then: 

di / / \fiy', a{y') - zn) - S,\^z% dzNdy' = 

J \y' —x'\<r ^max{0,ii(a;)— r} 



= (i4 

J I 

< Cd^r^ 



r+d(x) 

max{0, d(x)—r} \J\y'—x'\ 
r+d{x) 



zn I ! l/(y', fl(y') - zn) - i?dy' I dzM < 

\J\y'-x'\<r J 

dy' dz]\f < 



max{0, d{x)—r} \ l?^'^^;' |<r 



df ^ df a{y') 



dy' dyn dy' 



< Cdy / / \Vf\\y',a{y')-ZN)z% dy'dzN < C-^{l+Arr^ / 

Jmax{0,d{x)~r} J\y'-x'\<r "3 JB{x,r)nn 



(y)|V/|2a!y. 



In the above argument we made the following choices 



C = C{zn) ■= / f{y',a{y') - ZN)dy' r = ujj^\ / f{rz',a{rz') - ZN)dz' , 

\J\y'~x'\<r J J\z'-x'\<l 

and C being the Euclidean Poincare constant on the — 1 dimensional Euclidean ball of radius one. 

Since for any ^ e K, \ f - < 2\f - (.{zn)\'^ + 2\(,{zn) - in order to prove (fTTT|) in this case, it 
only remains to estimate the following term 

p pr+d{x) 

/ / \i{zN) - e?^^ dzNdy' = 

J \y' —x'\<r J uiax{0,d{x)—r} 



11 



dy'\ |C(^7v)-CT 



max{0, — r}°'~^^ 



\y'—x'\<r 
r+d{x) 



a + l 



r+d{x) 

max{0,d{x)—r} 



a + l 



max{0, d(x)~r} 



max{0, d{x) - r}"+^) dzN 



Thus, choosing ^ := ^(r + d{x)) above, we obtain by Holder inequahty 

fr+d(x) 



/ , / 

J \v'—x'\<r Jn 



y'—x'\<r J ma,x{0, d(x)—r} 



< 



a + l 



r+d{x) 



\y'—x'\<.r J max{0, d{x)—r} 



\^{zN)-^\^ztf dzNdy' j X 



r+d{x) 
\y'—x'\<r J max{0,d(a;)— r} 



dazN) 



dzjsf 



max{0,d(2;) — r}"^"*^ 



Since 

d^ZN) 



dzN 
hence 



-1 
N-1 



\z'-x'\<l 



Of 



dvN 



d^ZN) 



dz 



N 



) - ZN)dz' = 






df 


J\y'-x'\<r 





f 


df 


l\y'-x'\<r 





dz]sidy' 



{y',a{y') - ZN)dy' , 



{y',a{y') - ZN)dy' . 



Thus, since d{x) < -/r, we obtain ((TTT|) with constant Cp := 2^(1 + A)" (C + ^ 



We next prove the following local weighted Moser inequality: 



4(7+1)^ 



// 



Theorem 2.6 (Local weighted Moser inequality) Let a > 0, N > 2 and Q C be a smooth 
bounded domain. Then there exist positive constants Cm o-nd R = R{a, il) such that for any v > N + a, 
xefl,0<r<R and f G CQ°{B{x,r)) we have 



[ <i"(y)|/(y)|'(^+^)(iy<CMrV( 

Jl3{x,r)nQ 



x,r) V 



d^{y)\Vffdy 



B{x,r)r\^ 



B{a;,r)nQ 



Proof: Let us first consider the case where d[x) > ^r. By the standard Moser inequality, there exists a 
positive constant C such that for any x £ Q and any > if > 3 or any > 2 if = 2, the following 
holds true 



B{x,r) 



\fiy)\^('+^)dy<Cr' 



f V 



B(x,t) 



B(x,r) 



\f\^dy] , y fGC^{Bix,r)) 



(see for example Section 2.1.3 in |S('2j ). Thus we have 

2 

f d"(y)|/(y)|2(i+i)dy < {d{x) + rrCr\-"-^ ( [ \Vffdy] ( [ \f\'dy] ' < 

JB(x,r) \JB(x,r) / \J B{x,r) 
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2 



< Cr^ (tT^ ] (^""(dix) - rrri [ I riy)\Vf\'dy ] [ I riy)\ffdy 1 < 

\d{x)-rj \JB{x,r) J \JB{x,r) J 



< 



CMr'Vix,r)-"- ( [ d-iy)\Vf\'dy] ( [ d''iy)\f\'dy 

\JB{x,r) J \JB{x,r) 



2 



where Cm := C ^1 + j ^"^ ^^^'^ '^^ constant appearing in the volume estimate in Lemma 

[Qwhen A = 0. 

Let us now consider the case where d{x) < ^r. Then we claim the following local weighted Sobolev 
inequality: there exist positive constants Cs and R = R{a, 17) such that for any x € < r < i?, 
satisfying d{x) < ^yr, and any / G C^{B{x,r)), we have 



iV+ci-2 

2(iV+c) _ \ '^+" , 



d"(y)|/(y)|ivT^dy < C5 / d"(y)|V/|%. (2.12) 

B{x,r)nn J JB{x,r)nn 

If we accept (|2.12() . then the result follows, with Cm = Cs by means of Holder inequality, in fact we 
have 

/ d^{y)f^^^^^dy<{ f d"(y)/TOdy| ^ (/ d^{y)fdy\^ 

JB{x,r)nn \JB(x,r)nn J \JB{x,r)nn J 

as well as for any i' > N + a 

2(v-N-a) 

d"(y)/2(i+a)dy I ( / d''iy)fdy] <( f d^{y)f^^^^)dy\ y(x,r)^Wir 

B{x,r)r\Q. J \JB(x,r)nn J \JB(x,r)nn J 

In the sequel we will give the proof of (ITT^ . We will follow closely the argument of |KMT2| . If y c 
is any bounded domain and u £ C°°{V) then it is well known that 

Sn\\u\\ N < \ \Vu\\li(v) + lkl|Li(9y) , 

[V ) 

_ 1 

where := Niri ^T(l + y)] ^ (see p. 189 in [M])- Let us fix from now on that V := B{x,r) n ft, and 
let us apply the above inequality to u := d^f, for any / G C^{B{x,r)) and any a > 0. Thus we get 



SN\\d'=^f\\N < / (|V/|d" + aa!'^-i|Vd||/|)dy. 

Let us remark at this point that boundary terms on dfl are zero due to the presence of the weight d"^, 
a > 0. To estimate the last term of the right hand side, we will make use of an integration by parts, 
noting that Vd • Vd = 1 a.e.; that is we have: 

a [ d''-^\f\dy = a [ Vd -Vd d''-^\f\dy = [ Vd'' ■ Vd\f\dy = 
Jv Jv Jv 

= - [ d^Adlfldy- [ d''Vd-V\f\dy+ [ d^Vd ■ v \f\ dS. 

Jv Jv JdV 
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Under our smoothness assumption on 17 we have that \dAd\ < cq6 in fi^ for 5 smah, say < 6 < Sq, 
and for some positive constant cq independent of 6 {5q, cq depending on 0,). Now, if d(x) + r < 6, that is 
if r < , we have that V C ^Is and it follows that 

a [ d''-'\f\dy < coS [ d''-'\f\dy+ [ d''\Vf\dy, 
Jv Jv Jv 

[ d''-'\f\dy <{a- co5)-' [ d''\Vf\dy . (2.13) 
Jv Jv 

Consequently for any r G (0, -R(a, f])), R{a,Q,) := ;^min{(5o, ^} and any a > the following inequality 
is true 

a 



hence 



+ 1^ j^d-\Vf\dy . (2.14) 



LTT^(V) \a-CQ5Q 

To proceed we will use the following interpolation inequality (cf. Lemma 4.1 of |FMT2j ^ . 

Il>fll ^ ^(g- l) |l^a^ll , l) ||^a-l^|| 

\\d tWLoiV) < \\d /ll^^(^) + \\d , 



Vl<g<— b:=a-l + - -N , a>0. (2.15) 

^ - N -I q ^ ^ 

^From (|2.13j) and (|2.14j) . we get for any a,b,q as above the following inequality 

\\d'f\\LHV)<Ci\\d''Vf\\mv) , (2.16) 
where := f + lU ^^-^ 



q Sn \a-coSo J Q \a~co5o ^ 

Let US now apply inequality ()2.16|) to |/|'' instead of /, for s := ^ + 1, q := ^ , b := Bs. Due to 1)2. 15() 
we have a = 6 + 1 - = ^ + A, where A := B + I - ^^N. In this way we obtain 



<Ci{^ + l) [ld^^^\f\^\^f\dy^< 



1 , ,1 



<C2( (^j^d'^\vf\W 



where Cs := Ci (^f + lj. 

After simplifying we see that we have proved the following: there exists R = R{^^ + A, Vl) such that 
for all < r < -R and all x G O with d{x) < ^r, there holds 



/ <C [ d^^iy)\Vf\^dy , 

JB(x,r)r\n I JB(x,r)r\n 



iB{x,r)r\n J JB{x,r)r\n 

for any iV > 2 and any / e Co°(5(x, r)) under the following conditions A:= B + l- ^^N, ^ + ^4 > 0, 

2<Q<ooifiV = 2, 2<Q<^ifA^>3; here = = C^{N, Q, B, cq, 6o). 

Taking A = Q := and i? := ^ we deduce the local weighted Sobolev inequality ()2.12|1 with 

Cs = Cs{N, a, cq,Sq) and this completes the proof of Theorem 12.61 

// 
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Remark 2.7 Note that the upper bound for the length of the "balls" in the local weighted Moser inequality, 
denoted by R{a,0,), goes to zero as a tends to zero. 

Remark 2.8 Let us note that when N = 1, the corresponding analogue of the local weighted Sobolev 
inequality \2.1^) when = (— 1, 1) is the following one 

/ (1 - \y\T\f\'^{y)dy < Cs r" +— / (1 - \y\r\f'\\y)dy , 

for any f G Cq^{x — r, x + r), and any q>2if0<a<l and 2 < q < "^^^^^^ if a > 1. Consequently 
Theorem \1.5\ as well as its consequences can be also stated for N = 1. 

^From the results within this subsection, we will now deduce a new parabolic Harnack inequality up 
to the boundary for the doubly degenerate elliptic operator defined in H2.1() . To this end let us first 
make precise the notion of a weak solution 

Definition 2.9 By a solution v{y,t) to vt = —L^v in Q := {B{x,r) n 0} x (0, r^), we mean a function 
V e C^{{0,r'^);L'^{B{x,r)nn,\y\^ d''{y)dy)) nC^{{0,r^);H^{B{x,r) nn,\y\^ d'^{y)dy)) such that for any 
$ G C°((0, r^); Co°(i3(x, r) n ^l)) and any < ti < t2 < r"^ we have 

r / {\y\^d''iy)vt^ + \y\^d''{y)VvV^}dydt = . (2.17) 

Jt^ JB(x.r)nn 



Then we have 



Theorem 2.10 Let a > 1, A'' > 2, A G [2 — A'^, 0] and C R^ be a smooth bounded domain containing 
the origin. Then there exist positive constants Ch o-nd R = R{0,) such that for xGQ, 0<r<R and 
for any positive solution v{y,t) of ^ = j;^jp^^^div{\y\'^d°'{y)'Vv) in {B{x,r) n 17} x (0, r^), the following 
estimate holds true 



ess 



'''Piy,t)<,{B{xr)nn}>,{rl,i^f^y^i)<CH ess ^nf ^y^^^^^^^^^r^^^^^^s^,^^,^v{y, 



In order to prove the parabolic Harnack inequality in Theorem 12.101 we use the Moser iteration 
technique as adapted to degenerate elliptic operators in |FKSj . |CSj as well as |GSCj . In this approach 
one inserts in the weak form of the equation vt = —L^v suitable test functions $. One of the key ideas 
is to use test functions $ of the form rj'^v'^, where v is the weak solution of the equation, is a cut off 
function and g € R. To this end one has to check that q'^v'^ is in the right space of test function. In this 
direction the following density theorem is crucial. 

Tlieorem 2.11 Let N > 2 and C R^ be a smooth bounded domain. Then for any a > 1 

H\n, d'^iy) dy)=H^{n,d^iy) dy) . 

In particular for any a > 1, the set C^(0) is dense in H^{VL,d°'{y)dy). 

Here H^{Vl,d°'{y)dy) denotes the set {v = v{y) : j^d°'{y){v'^ + \Vv\'^)dy < oo}, the corresponding 
norm being defined in 1)1.12(1 . 

We are now ready to prove the density theorem. 

Proof Let us prove here the result when a = 1. We refer to Proposition 9.10 in "K" for the case a > 1, 
even though our proof with some minor changes, can also cover this range. 
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First of all from Theorem 7.2 in ^ it is known that the set C°°(n) is dense in H^{Q, d{y) dy). Thus 
for any v G H^{i^,d{y) dy) there exists Vm G C°°(r2) such that for any e > we have \\v — Vm\\H^ < e if 
m > m(e). Let us choose w := fm(e) and let us define, for A: > 1, the following function 



if d{x) < 4 



1 + ^ ^f^<dix)<l 



In(fc) 'J W ^ ""^-^^ ^ k ' 
if d{x) > i . 



Then Wk ■= Wipk G Cq'^{^), moreover we have 



\W - Wk\ 



|„i = < 2 / {w^ + \Vw\^){l-ipk?d{y) dy + 2 [ w^\Vvk\^d{y) dy < 

JQ Jn 

<2 [ {w^ + \Vw\^)d{y)dy + 2 [ dy . 

JdivXl Jj^<d(v)<^ d(y)(\n(k]y 



IdivXi"'^ ' ~J^<d{y)<j: d{y){ln{k)f 

Now as A: ^ oo the right hand side goes to zero, this proves the Theorem. 

// 

The above Theorem allows us to take the cut off function r] in C^{B{x,r)) instead of taking it as 
usual in Cq°{B{x, r) n 17). Clearly the two function spaces differ only if the "ball" intersects the boundary 
of Q. To explain what are the appropriate modifications of the standard iteration argument by Moser, we 
now present in detail the first step, which is the mean value inequality for any positive local subsolution 
of the equation vt = —L^v. 

Theorem 2.12 Let a>l, N >2, AG [2 — A'^, 0] and C be a smooth bounded domain containing 
the origin. Then there exist positive constants C and R{0,) such that for x £ 0,, < r < R{0,) and for 
any positive subsolution v{y,t) of vt — \y\\d"{y) div{\y\^d'^{y)Vv) = in {B{x,r) n 0} x (0, r^) we have the 
estimate 

C f 

ess sup, ^, r„, r^ ^i / r2 ox v"^ (v , t) < " / Ivl^ d"' (y) v'^(y,t)dydt . 

Proof: We will only prove the result in the non standard case in which the "ball" B{x, r) intersects 
the boundary of 17; we refer to |MT2j as well as to |GSCj for details in the other case. Similarly to 
Definition 12 . 91 we define a subsolution v{y,t) to be a function in C^((0, r^); r) n J7, \y\^d°'{y) dy))n 

C^{{0,r'^);H^{B{x,r) n 17, dy)) such that the following holds true 

/ / + |y|V(y)VwV$}dydt < 0, <^> e C\{0,r^);C^{B{x,r) nfl)) > 0. 

Hence in particular we have also 



[ {\y\^d'^{y)vt^ + \y\^d''{y)VvV^}dy < 0, V $ G C^{B{x, r) n 17) , $ > 0. 

JB{x,r)r\n 



g-l 



Let us define for any q,M > 1 the following functions G{z) = z'^ \i z < M and G{z) = M'^ + q{z — M)M 
if z > M and H{z) > by H'{z) = y^G^, H{0) = 0; note that G{z) < zG'{z) as weh as H{z) < 
zH'{z). Due to Theorem 12 . 1 II there exists a sequence of functions Vm in C°°{B{x,r) n 17) having compact 
support in 17 such that Vm — > f in H^{B{x,r) n 17,d"(y) dy) as m ^ +00; whence due to ()2.4|) also in 
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H^{B{x,r) n Q, dy). Hence for any rj G C^{B{x,r)) and m > I the function $ := ■rfG{vm) is 

an admissible test function, that is the following holds true 

{\y\''dP{yWG{vm)vt + \y\^(P{y)VvV{ifG{v^))]dy < . 

B{x,r)nn 

Passing to the limit as m ^ +00 we get 

{\y\''d''{yWG{v)vt + \y\^d''iy)VvViv''G{v))}dy < , V r? G C7o^(e(x,r)) . 



B{x,r)nn 

This is the standard starting point in Moser iteration technique apart from the fact that the cut off 
function r] does not be necessarily zero on di}, this is crucial. Then by Schwarz inequality we get 

/ {\y\^d''{y)v^G{v)vt + \y\^d''{y)\Vv\^G'{v)ij''}dy<C I \y\>^d^{y)\V7^\\^G'{v)dy 
JBix,r)nn JB{x,r)nn 

thus also that 

{\y\^d"iyWG{v)vt + |y|^d-(y)|V(77i/(t;))p}dy <C [ \y\^d''{y)\Vv\^v^G'iv)dy . 

B{x,r)r\D, JB{x,r)nn 

For any smooth function x of the time variable t, we easily get 

1 / \y\^d''{y){vxF{v)fdy + [ \y\^d''{y)\V{vH{v))\'dy < 

JB{x,r)nn JB{x,r)nn 

<Cx(xl|Vr/||ioo(R„) + ||x'||L^(R)) / \y\^d''{yyG'{v)dy; 

here F{z) is such that 2F{z)F'{z) = G{z). For i < s < s' < 1 we choose as usual x such that < x ^ 1; 
X = in (— cxD, r^(l — s')), x = 1 in (r^(l — s), 00), moreover if ^ G C(j"(0, 1) be a nonnegative non increasing 
function such that ^(z) = 1 if z < s and ^(z) = if z > s', we define, making use of local coordinates, the 
following cut off function r/(y) := ^ (y^^—f^^ S, ^ l"*-^ )~yN-dix)\ \ xhen clearly ||Vr7||ioo(j^n') < ^^^^^.^ and 



IIx'IIl°=(r) < 

Integrating our inequality over (0,t), with t G (r2(l - s),r'^) we obtain 

sup/ \y\^d"{y){vF{v)fdy+ [ \y\^d''iy)\V{r^H{v))\^dydt < 

t€J JB(x,r)nn J{B(x,r)nn}x{r2(l-s),r2) 

^ 2^ ^2 / \y\^d^iyyG'iv)dydt . 

r [S -S) J{B{x,s'r)nn}x{r^{l-s'),r2) 

Making once again use of Theorem 12.111 we note that we can apply the local weighted Moser inequality 
in Theorem 12.61 to the function / := r]H(v) thus obtaining 

/ \y\^d''{y){r]H{v)f(^^^hydt < 

J{B{x,r)nn}x{r^{l-s),r^) 



< ^7^^^ ( / \y\^d''{y)v^G'{v)dydt 

'{B(x,s'r)nf7}x(r2(l-s'),r2) 



N + a 
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Let us now denote by 7 := 1 + thus as M tends to infinity we have for p := g + 1 



/ \y\^d''{y)v^dydt < f ( / \y\^d''{y)vPdydA 

J{B{x,sr)nn}x(r^(l~s),r^) r [S - S) yj {t3(x,s'r)nn}x(r^{l-s'),r'^) 

Thus due to Lemma 12.21 also that 

V{x,sr)-\r^s)-^ [ \y\^d'^{y)v^ dydt < 

J{B(a;,sr)nf7}x(r2(l-s),r2) 



<c(^^\ (v{x,s'ry\rh')-^ [ \y\^ d'' {y)vP dydt 

\S - S J y J {B{x,s'r)r\Q.}x(T'^(l~s'),r^) ^ 

Take now p = pi := 27*, s = ^j+i and s' = Oi where 6i := „?^^^,^ then if we denote by : 



1 



2 / ^ \ 

{B(x,s'r)nn}x(r2(l-s'),r2) 
2{i+l) 

y(x,6'ir)~^(r26'j)"^ /{B(a-,eir)nn}x(r2(i-e,),r2) \y\^d"'iy)v^'dydt'j the above inequahty can be restated as 
follows /(« + 1) < C{i)I{i). Thus since one can show that the product of C{i) for all i > is finite, 
we obtain /(oo) < {HSo ^(01 this completes the proof of the proposition. To this end the choice 
:= min{/3, 0)} can be made, here (3 and fi) are the constants appearing respectively in the 
local representation of 90 and in Theorem 12.61 when a := 1. 

// 

Theorem 12.61 corresponds to the local weighted Moser inequality needed in the proof of the parabolic 
Harnack inequality up to the boundary stated in Theorem 11.51 The local weighted Moser inequality 
involved in the proof of Theorem 12.101 differs from Theorem 12.61 onlv if d{x) > 'jr, N > 3, X ^ 0, and in 
this case it reads as follows 

Theorem 2.13 Let N > 3, A G [2— A^, 0) and ft C be a smooth bounded domain containing the origin. 
Then there exist a positive constant Cm such that for any i^>N,xGfl, r>0 and f E CQ^{B{x,r)) we 
have 

2 

[ \y\'\f{y)f^'^^^dy < CMr2(r^(|x| + r)')-'' ( [ \y\^\Vf\'dy] ( [ \y\Wdy] " • 

JB{x,r) \JB(x,r) ) \JB(x,r) ) 

Proof : By Holder inequality the result easily follows with Cm '■= Cs as soon as the following local 
weighted Sobolev inequality holds true 

iV-2 

> 2JV \ ^ 2|A| r \ n 

y\^\f{y)\~'dy\ <Cs{\x\+r)— \y\^\Vf\^dy (2.18) 

B(x,r) J JB{x,t) 

(we refer to the proof of Theorem 12.61 where a similar argument is used). Let us first prove the above 
inequality for any AG (2 — A^, 0). As a consequence of the Caffarelli Kohn Nirenberg inequality (e.g. see 
Corollary 2 in Section 2.1.6 of |M1)) the following holds true 

N-2 

If f^2\y\w^2dy] <C [ \Vf\^\y\^dy, yfeC^{B{x,r)) 

\JB{x,r) J JB{x,r) 

and for some positive constant C independent of x and r. Whence also that 



7. 



N^2 2|A| 
iV / \ iV 



f^^\y\^dy\ <C[ sup / iV/Plyj^dy < C7(|x| + r)^ / \V f\''\y\^dy . 

B{x,r) I \yeB{x,r) I J B{x,r) J B{x,r) 
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Let us now prove the result for A = 2 — A^. To this end let us apply Proposition 13. II to 0, = B(0, 1) with 

1 

D = . Then there exists a positive constant C such that 

N-2 

J 3(0,1) \J 3(0,1) \DJ J 

here X{t) = j^^, t £ (0, 1]. Now let us take v{x) := / (-1) for any / £ Cq°{B{0, R)) then from above we 
have 



JV-2 
JV 



Jb(o,r) \Jb(o,r) \^^J ) 

Then if y G B{x,r) clearly y G B{0, \x\ +r), thus if we take R := \x\ +r and / G C^{B{x,r)) from above 
we have 

JV-2 

/ wrnvr-dy >c([ /^ir^x^ ( m A ^ > 

JB(x,r) \JB(x,r) \^ ti ) \ 



2N 



N-2 

N 



iB(x,r) J \yeB(x,r) \DR 

Whence the claim easily follows as soon as we prove that 



\y&B(x,r) 

This is indeed the case in fact we have 



2{N-2) 

sup \y\X-^{^]] <Cs{\x\+r)'-^ 



JV-l 



DR 



N-l 

sup (^) < sup \y\ (l - In (^^^ 

yeB(x,r) \^^J 0<|s/|<|x|+r V \Uri, 



N-l 



(thus using the fact that the function ip{t) = t (l — In (;^)) is an increasing function for t G [0, R] if 
D and R are as above) 

= {\x\ +r) 1 - In ' ' = +r)(l + ln(Z)))~ = (|x| +r) 



DR J J M I < /V ■ V // Ml ^ \N -2 

This completes the proof of Theorem 12.131 

// 

To state the heat kernel estimates following from Theorem 12.101 we introduce some notation. The 
operator is defined for a > 1 and A G [2 — A^, 0] in L'^{Q,\x\'^d°'{x) dx) as the generator of the 
symmetric form 

'^a[vi^V2] := / \x\^d°'{x)VviVv2 dx , 
Jn 



namely 



D{L^) := G H^{n, dx) : --^^^p^^^div{\x\^d''{x)Vv) G L\n, |x|^d"(x) dx) } , 

L^v := -^^^^^^^div{\x\^d^{x)Vv) for any v G D(L^) , 
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where Hq{Q, dx) denotes the closure of (7^(17) in the norm 

\v\\hi_ ■■= \x\^d'^{x){\Vv\'^ + v'^) dxY . (2.19) 



Then is a nonnegative self-adjoint operator on L^(J7, \y\^d'^{y)dy) such that for every t > 0, e 
has a integral kernel, that is e~^°'^vo{x) := J^l^{t, x,y)vo{y)\y\^d°'{y)dy; here l^{t,x,y) is called the heat 
kernel of L^. The existence of l^{t, x, y) can be proved arguing as in |DSlj : that is, using a global Sobolev 



inequality on 0, which can be easily deduced from its local version (|2.12|) as well as (|2.18|) . by means of 
the partition of unity as in |K] . 

Then from the parabolic Harnack inequality in Theorem l2.ini the following sharp two-sided heat kernel 
estimate can be easily deduced: 

Theorem 2.14 Let a > 1, N >2, AG [2 — A'^, 0] and C R'^ be a smooth bounded domain containing 
the origin. Then there exist positive constants Ci, C2, with Ci < C2, and T > depending on $7 such that 

[t2 d2{x)d2{y) J 

< mm (4, M±4)l(^±^| e-^^^ 
\t2 d2(x)d2{y) J 

for all x,y £ il. and < t < T . 

Proof of Theorem \2.14\ ' Using the mean value estimate for subsolutions as in Theorem 12.121 and the 
parabolic Harnack inequality of Theorem I2.1UI and arguing as in Theorems 5.2.10, 5.4.10 and 5.4.11 in 
!2] we are lead to the following Li-Yau type estimate 



-<l^{t,x,y)< 



V{x,Vi)^V{y,Vif2 " ' ' V{x,Vi)^V{y,Vi) 



1 ) 

2 



for all x,y G and < i < T; where Ci, C2 are two positive constants with Ci < C2, and T > depends 
on n. From this the result follows using the volume estimate in Lemma 12.21 

// 

Using all the machinery we have produced in this section we can handle more general operators than 
the one in Theorems 12.101 and 12.141 Thus, consider the operator 

1 ^ a / . d 



■= - MA;^ E 7f;r. a,,ix)\x\'d-{x)— , (2.20) 



\x\'^d°'(x) dxi \ dxi 

1,3=1 ■' 

where (oj,j(2;))„xn ^ measurable symmetric uniformly elliptic matrix. The operator L\ is defined for 
a > 1 and AG [2 — A^, 0] in L^(r2, |a;|'^(i"(2;) dx) as the generator of the symmetric form 

N 



^i[vi,V2]:=Y^ / |x|^d"(x)aij(x) — — . 



Then the existence of a heat kernel l^{t,x,y) follows as in |DSlj . and we have 
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Theorem 2.15 Let a > 1, N >2, AG [2 — A^, 0] and 17 C be a smooth bounded domain containing 
the origin. Then there exist positive constants Ci, C2, with Ci < C2, and T > depending on Vt such that 

^ (|x| + ^/^)^(|^/| + Vt)^ ) c,i^ 

t2 d-2{x)d2{y) J 

/or a// x,y € Q and < t < T . 

Remark 2.16 A parabolic Harnack inequality up to the boundary similar to the one of Theorem A2.1(]\ can 
be stated under the same assumptions of Theorem \2.15\ for the more general operator L^. 

3 Critical point singularity 

In this section we establish a new Improved Hardy inequahty (Theorem I3.2|) and then we give the proofs 
of Theorem II . 1 1 and Theorem 11.21 The structure of this section is as fohows. 

In Subsection 3.1 we first deduce the improved Hardy inequahty and then the global in time pointwise 
upper bound for the heal kernel of the Schrodinger operator —A — ((A^ — 2)^/4)|x|~^, which is sharp when 
X and y are close to the boundary (see Theorem 13. 4() : then, due to an argument contained in |Dlj . we 
complete the proof of Theorem 11.21 proving the sharp lower bound for time large enough. 

The proof of Theorem ll.ll is finally completed in Subsection 3.2, using the parabolic Harnack inequality 
up to the boundary of Theorem l2.ini 

3.1 Boundary upper bounds and complete sharp description of the heat kernel for 
large values of time 

We first recall the following improved Hardy-Sobolev inequality stated in Theorem A in |FTj (see also 
inequahty (3.3) in [BFT2j) 

Proposition 3.1 For N > 3, let Q C R^ be a smooth bounded domain containing the origin and D > 
sup^gQ Then there exists a positive constant C such that 

I \Vv\'^\x\'^'^dx>c([ v^\x\-^X^^ 

for any v G C^{^); here X{t) = t G (0, 1]. 

We next state a new result, the proof of which will be given later on. 

Theorem 3.2 (Improved Hardy inequality) Let C R^, N > 3, be a smooth bounded domain 
containing the origin. Then there exists a constant C = C{^) G (0, |] such that 

The positive constant C(i7) can be taken to be exactly ^ for all domains satisfying the following condition 

-div{\x\^~^Vd{x))>^ a.e.inn. (3.2) 



Ci min 



< C2 min ■ 
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For example when Q = B{0,R), for arbitrary R > 0, condition H3.2jl is satisfied. Consequently, in this 
case the Hardy inequality involving the Schrodinger operator having critical singularity at the origin can 
be improved exactly by the inverse-square potential having critical singularity at the boundary. 

As a consequence of Proposition 13.11 and of the improved Hardy inequality of Theorem 13.21 the 
following logarithmic Sobolev inequality can be easily obtained: 

Theorem 3.3 (Logarithmic Hardy Sobolev inequality) For N > 3, let Q, C be a smooth 
bounded domain containing the origin. Then for any u G Co°(r2 \ {0}), n > 0, and any e> Q we have 

f uHog{ \dx<ej {\Vu? -^-^^-^u')dx+(K,-^\oge)\\u\\l- (3.3) 

\\\u\\2\x\ — d{x) ) Jn\ 4|a;|^ J \ 4: J 

here is a positive constant independent of e and \ \u\\2 ■= [J^u^dx)^ . 



Then using Gross theorem of logarithmic Sobolev inequalities, as adapted by Davies and Simon (see 
Theorem 2.2.7 in we will show the following global in time pointwise upper bound for the heat 

kernel: 

Theorem 3.4 For N > 3, let Q C be a smooth bounded domain containing the origin. Then there 
exists a positive constant C such that 

k(t,x,y) < C^^^My) e-^^*, yx,yen,t>0. 

Let us first prove the logarithmic Hardy Sobolev inequality 1)3. 3() . 

Proof of Theorem \cl.c\ As a first step we claim that the following logarithmic Hardy Sobolev inequality 
holds true: 

J u^{- log d{x)) dx<ej Vnj^ - dx + (-^^i - ^ log | \u\ \l , (3.4) 

for any u E Cq^{^), n > 0, and any e > 0; here Ki is a positive constant independent of e. 

To see this let us first suppose that the nonnegative function u G Cg°(r2) is such that ||m||2 = 1- We 
then have 



/ u^{-\ogd{x)) dx = ]- [ u'^{\ogd{x) ^) dx <^log( [ — 
Jn ^ Jn 2 Vic d{x^ 

< i log f / ( \Vu\' - dx 



u^dx I < 



2 7nV 4|x|2 

here we have used first Jensen's inequality and then the improved Hardy inequality (|3.1|) . For a general 
nonnegative u G Cq^{^) we apply the above inequality to the function tt^^TT' ^° 



/ 



u^{-\ogd{x)) dx < ^||u|||log (j^iil y {\^'^? - 



Since logz < z for any z > 0, then also logy < e2Cy — log (e2C), for any e > 0; whence from this we 
deduce with Ki := \ log(2^). 



22 



We will next show the following logarithmic Hardy Sobolev inequality: 



^'logf ,, ,, ,2-N ]dx<e j ("iVnl^ - .f ' ^') dx + - jloge ) ||n||i , (3.5) 



U 2 F 2 



for any u G C^(i7 \ {0}), n > 0, and any e > 0; here K2 is a positive constant independent of e. 
By Proposition 13. II it follows easily that there exists a positive constant C such that 



JV-2 

\Vv\'^\x\'^-^dx>C[^j v^\x\'^-^dx^ , (3.6) 

for any v E C^{Q) (this is inequality (4.12) in [BFT2I). Whence we claim that the following logarithmic 
Sobolev inequality holds true: 



J y2^og(^-^^ |x|2-^dx<ey |Vt;|Vl'"'^f^^+ (^2- jloge) 



12 ! 



(3.7) 



for any v G C^{^}), v > 0, and any e > 0; here K2 is a positive constant independent of e and ||f||2 : = 

(Jj^ f ^|xp^^(ix) ^ . To see this let us first suppose that the nonnegative function v G C^{^) is such that 
||f||2 = 1. We then have 

t)^log(f) \x\'^~^ dx = ^ ^ ^ j v^\og{v^^^ |x|^^^ dx < ^ ^ ^ log (^j v^^^"^ \x\^~^ dx 



N-2 



= ^log(^j^^il^ |x|2-^dx) " <^log(c-i j^|V^;|2|x|2-^dx) ; 

here we have used first Jensen's inequality and then the improved Hardy-Sobolev inequality (|3.6|) . For a 
general nonnegative v G CQ^{Vi) we apply the above inequality to the function ypjjjj^, to get 

y^log(^-^^ \x\''-'' dx<^\\v\\l\og(^^J |Vz;p|xp-^dx 

Since log 2; < z for any 2; > 0, then also logy < e^j^y — log (e-jy)' ^'^'^ ™y e > 0; whence from this we 
deduce (jSH with K2 := f log (^). 

Inequality 1)3. 7p implies (|3.5p via the following change of variables u := t;|x|^~. Finally from 1)3.4(1 
and (|3.5|) . the logarithmic Hardy Sobolev inequality i|3.3(l easily follows with constant K^, := Ki + i^2 + 
^log2. 

// 

We are now ready to give the proof of Theorem 13.41 

Proof of Theorem \3.4\ Let us define, as in Section 2 of |D2], the operator K := U'^iK - Xi)U, 
U : L'^{Q, ip\dx) —>■ L'^{Q) being the unitary operator Uw := ipiw, thus K := —-^div{^p\V). Here 991 > 

denotes the first eigenfunction and Ai > the first eigenvalue corresponding to the Dirichlet problem 
— A931 — ifi = Xi^pi in 0, (^1 = on dO., normalized in such a way that ^'lix) dx = 1. Due to 



the results in Lemma 7 in |DDj and using Theorem 7.1 in |DSlj on one hand and elliptic regularity on 
the other, there exist two positive constants ci,C2 such that 

2 — iV 2-N 

ci\x\ 2 d{x) < Lpi{x) < C2\x\ 2 d{x), V X G $7 . (3.8) 
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^Prom this and (|3.3j) we deduce the following logarithmic Sobolev inequality 

X ( IRI2 ) '^1^^ - ^ ^ ^^'^ + (-^4- ^^^ loge^ 1 1 will , (3.9) 

for any w G Cq°{Q, \ {0}), w > 0, and any e > 0; where K4 := — logci and ||w||2 := (/q w'^ipf dx)^ . 
Let us remark that only the lower bound in estimate (|3.8|1 was used. 

^From now on one can use the standard approach of |D4j to complete the proof of the theorem. Here 
are some details for the convenience of the reader. 

As a first step we claim that the following logarithmic Sobolev inequalities holds true: 



If^w^^og ^ldx<e?^< Kw,wP-' >LW, d.) + {ka - ^^loge) ||^z;||^ (3.10) 

for any w G C^(J7 \ {0}), w > 0, and any e > 0, p > 2. To see this we apply inequality (|3.9() to w^; 
whence due to the fact that 

[■ 2 /• 2 

\yw'i\W dx = ^ J^wP-^\VwWl dx = j^^—^ < Vw,VwP-^ >L^{n,vl dx)< f < Kw^w^-^ >L^(n,vl dx) > 

since 2{p-i) — 1 if P — 2; the claim follows. 

Let Hq(Q, ifl dx) be the closure of Cg°(r2) with respect to the norm 

\\w\\h^ '■= \ ilVwl"^ cpl + w'^ipj) dx 

as one can easily prove this is also the closure of C^{^} \ {0}) with respect to the same norm. Then 
to the operator K defined in the domain D{K) = {w G Hq{Q.,(p\ dx) : Kw G L?'{Q.,ip\ dx)} it is 
naturally associated the bilinear symmetric form defined as follows K,[wi,W2\ ■=< K'Wi,W2 >L^(n^,pj dx)~ 
Vwi'S/w2 dx, which is a Dirichlet form. Whence Lemma 1.3.4 and Theorems 1.3.2 and 1.3.3 in |D4j 
implies that e~^*, which is an analytic contraction semigroup in L'^{Q, tp\ dx), is also positivity preserving 
and a contraction semigroup in U\Q, Lp\ dx) for any 1 < p < 00. As a consequence for any t > and any 
p>2 

e-^*[L2(0, ^\ dx) n L°^{Vt)]+ C [Hl{^, ifj dx) n L^i^, ipj dx) n L°°(rj)]+ ; 
where we denote by [E]'^ the subset of positive functions in the space E. 

Thus by density argument the logarithmic Sobolev inequality (|3.1U|) . more generally applies to any 
function in Ut>o e-^^[L^{n,ipldx) H L^{Q)]+ . 

This means that Theorem 2.2.7 in |D4j can be applied, in the same way as in Corollary 2.2.8 in |D4j. to 
the operator K; whence obtaining that 

||e-^*||2^oo <Ct-— , 

and by duality that 

||e-^*||i-^2 <Ct-— , 

that is 

lle-^'^lli^oo < Ct~— . 
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Here we use the following notation: 



-Kt\ 

0<ll/ll9<l 



\e-^'f{x)\\, 



" sup ii^^^^ii^ 



where ||/||q := (/^^ l/l'^'/'i dx)'' . This implies, by Dunford-Pettis theorem, that the semigroup e is 
indeed a semigroup of integral operators; that is a heat kernel k{t, x, y) associated to the semigroup e~^^ 
is well defined and satisfies the following pointwise upper bound k{t,x,y) < Cjt ^, for any x,y ^ Q 
and any t > 0. Theorem 13.41 then follows, due to the upper bound in (|3.8() and to the fact that, as a 
consequence of the unitary operator U, the heat kernels k{t, x, y) and k{t, x, y), corresponding respectively 
to K and K, satisfy the following equivalence 

k{t,x,y) = Mx)My) Ht,x,y)e-^^' . (3.11) 

// 

Remark 3.5 Applying Davies's method of exponential perturbation to the operator K (see Section 2 in 

\x — y\^ 

ID3\1 for details), the upper bound in Theorem\3.4\ can be improved by adding a factor c^e . 



Let us now deduce from the upper bound in Theorem 13.41 an analogous lower bound for time large 
enough, thus completing the proof of Theorem ll.2l We argue as in Theorem 6 of |Dlj (see also Proposition 
4 of |D2j ). we give the details here for the convenience of the reader. 

Proof of Theorem Making use of the same notation as in the proof of Theorem 13.41 the lower 
bound we want to prove corresponds to the statement k{t, x,y) > C for any x,y £ Q if t is large enough, 
C being some positive constant. 

For any / G L^{Q,ip1 dx), we clearly have 

f=<f,l>l + g, 

where < /, 1 >:=< /, 1 >L2(n,<fil dx)i ^-^d < g,l >= 0, since ipl{x)dx = 1. Thus, making use of the fact 
that by definition i^'l = we have 

e-^^V =< /, 1 > 1 + e-^*5, 

that is the semigroup e~^^f := e~^^f— < /, 1 > 1, to whom it is clearly associated the heat kernel 
k{t, X, y) — 1, is such that for any / G L^{Q., ip\ dx) 

where g = g{f) is a function in L^{Q, (p\dx) such that < 5, 1 >= 0. Thus, due to Theorem 13.41 

||e lli^oo S ||e lli^oo S C/t 2 J 
here C is some positive constant; this is equivalent to say that 

\k{t,x,y)-l\ <Ct-— , 
from which the claim easily follows for t large enough. 

// 

In the sequel we will give the proof of Theorem l3.21 We will use the following lemma whose proof will 
be postponed until the end of this subsection. 
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Lemma 3.6 For N > 3, let Q C be a smooth bounded domain containing the origin. Then there 
exists 60 > 0, such that 

for all < d < 6q; here il^ := {x G : dist(x, dQ) < 6}. 

We are now ready to prove the improved Hardy inequality. 

Proof of Theorem \3.}^ (i) Let us first prove the claim on any domain satisfying condition ^i.2\ . 
To this end let us define for any u G C^{Vl) as a new variable w := |x| 2 d~2(x)ti, obviously w G Hq{^). 
By direct computations we have 

2 — N ,_iV_-| 1 2-N 1 _1 2-iV 1 

Vii = — ^ — |x| 2 X (i2 u; + 2 -d 2\/d w + \x\ 2 d^Vw 

thus 

+ y ^- ^ ~ -w^ X Vd - (iV - 2)|x|"^ dw xVw + VdVw w dx . 

Vu;|2 d - X Vd - i^^^|x|-^ dxVw^ + -Vd Vw^ \x\^-A dx = 

III 2'' 2'' 2 y 

Vw;pd|x|2-^ - ^^|x|-^ ^ vd+ ^^^^diz;(|x|-^d x)^;^ _ ^rf^^(|^|2-iVy^)^2 j ^ 

= y (iVu-l^ d |x|2-^ - ^dit;(|x|2-^Vd)u;2^ dx > 0, 

due to condition H3.2|) on 0. Thus inequality ()3.1|) is proved with constant C(0) = ^ in any domain Q. 
satisfying condition (|3.2p . 

(ii) Let us prove indirectly the claim in the remaining case. To this end let us denote by H^ip., |xp~^dx) 
the closure of CQ°(r2) in the norm 



Whence 



-'2-JV 



HI „:= W (|V/|2 + /2)|x|2-^dx^ . (3.12) 



JV — ^ 

The improved Hardy inequality ()3.1() we are going to prove, in the new variable v := |x|^~u reads as 
follows 

|Vt;| Vl^-^dx >C I |x|2-^4^^ • 
Jo. 

Let us suppose that the improved Hardy inequality (|3.1() is false; whence let us suppose that the following 
holds true 

inf / |xp-^|V?jpdx = ; 
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thus there exists a sequence {vj}j>o in Hq{Q, ^dx) such that dx = 1, and 

\x\^~^ \SIvj\^dx ^ as j — > oo . (3.13) 
For any arbitrary function ip £ C^{0,), such that 99 = 1 in a neighborhood of the origin, we also have 

/ \x\^-''\V{ipvj)\^dx<2 [ |x|2-^(|Vt;,|V^|Vv;|27;2) dx < 
Jn Jn 

<C [ {\^Vj\^ + v]) dx<C [ \x\'^~^\Vvj\'^dx ^ as j ^ 00 . (3.14) 

Jn Jn 

Here we use the fact that the foUowing inequahty holds true 

[ \x\^~^fdx<C [ \x\'^-^\Vf\^dx, y f e H^{n,\x\^~^dx) . (3.15) 
Jn Jn 

Inequality ()3.15|) for example follows easily from inequality p.6() by Holder inequality. From estimate 
()3.14|) and inequality 1)3. 15(1 (applied to / := ipvj) we easily deduce that 

\2-N 2 2 n 
X\ if Vj — > 0, as J — > CXD, 

n 

or similarly (due to the fact that ip has compact support inside Q) that 

\x\^~^ip^^dx asj^cx). (3.16) 



We then compute 

1,2 



1 = f \x\^-^^dx = [ 
Jn Jn 



|^|2-jv (V^^i + (l-<^)^j)^ ^ 
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+ \x\^'''^{l-v)''-idx+ t \x\^~^{l-^f'ldx. 
d^ Jn d^ Jq d^ 

We observe that the first two terms in the last line tend to zero as j tends to infinity and therefore we 
obtain that 

\x\^'^il-<pf4^dx = l + o{l), asj^oo. (3.17) 

On the other hand we have that 

[ \x\'^-^\V[{l-^)vj]\^dx<2 [ \x\^~^\Vvj\'^dx + 2 [ \xf-^\V{ipvj)\^dx , 
Jn Jn Jn 

both terms in the right hand side going to zero as j tends to infinity due to (|3.13|) and (|3.14|) : whence we 
deduce that 

\x\^-^\V[{l- ip)vj]\^dx ^0, asj^oo. (3.18) 

Jn 

Since for any j > the function / := (1 — ip)vj is an element of Hq{Qs) for a suitable choice of the 
function ip (take it identically one in a subset containing ^1 \ il-s), by means of (|3.17|) and (|3.18|) we reach 
a contradiction with Lemma 13.61 thus proving the improved Hardy inequality. 
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// 

A similar improved Hardy inequality for a potential behaving like {{N — 2)^/4)|a;| ^ near the origin 
and exactly like (l/4)d~^(x) near the boundary is also shown without any geometric assumption on the 
domain Q (see Theorem 13. lUI below) . 

We next prove Lemma 13.61 One can consider it as a consequence of the following more general result. 

Lemma 3.7 For N > (i, let ^ C be a smooth bounded domain. Then there exists a positive constant 

N 

^0 = ^o(^); such that for any V G L~{Qsq) cif^d < 5 < (^o? have the following estimate 

j ('^"'^ ~ i^^) ^^-^j ^u^dx , Vug C^i^s) ; 
here c = c{5) — > oo as 5 — > and il^ := {x G $7 : dist(x, dO.) < 5}. 

Proof of Lemma Vj.fA: Let us choose V{x) := in Lemma 13.71 above and let us choose 5 small 

enough such that c{5) > 1 and < 6 < 5o, thus we have 



-^uA dx > / i^)' 



l^^r - Tm'^ ]dx> I ' u-'dx, (3.19) 

Us 



for any u £ C^{^ls). For any / G C^{i^s), u ■= f\x\ 2 will be in C^{^ls), moreover by easy computa- 
tions we have 



thus ()3.19|) can be restated as follows 



this proves the claim. 



// 



Whence it only remains to prove Lemma 13.71 Before doing so let us observe that inequality ()3.19|) . 
simply says that the improved Hardy inequality (|3.1|) indeed holds true with constant C{Q) = | whenever 
the support of the functions considered is contained in a neighborhood of the boundary. 

The proof of Lemma 13.71 makes use of the following improved Hardy-Sobolev inequality near the 
boundary stated in Theorem 3 of |FMTlj . we recall it here for the convenience of the reader: 

Proposition 3.8 For N > 3, let Q C be a smooth bounded domain. Then there exist positive 
constants 60 = 5o{0,) and C = C{N), such that 

N-2 

and any Q < 5 < 6q; here 0^ := {x G : dist(x, dVL) < 5}. 
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Let us focus here on the fact that in ProDosition l3.8l no convexity assumption on the domain 17 is made; 
this is due to the fact that we only consider functions whose supports are contained in a neighborhood of 
the boundary. 

Proof of the Lemma 7| " By Holder inequality we have 

_2_ N-2 

Vu^dx< f / V^dx] i / n~dx| < 



ns / \Jns 



the last step being due to Proposition 13.81 This proves the claim with constant 

c{6) := (^j V^dx^ ^ C{N) , 
which tends to infinity as 6 tends to zero due to the integrability assumption on V. 

// 

With some minor changes in the proof of Theorem 13.21 one can indeed prove the following improved 
Hardy inequality, which does not a priori requires the bounded domain Q to be smooth. 

Theorem 3.9 For N > 3, let Q C be a bounded domain containing the origin such that 

,2 



[ \Vu\'^dx >c[ ^dx, Vug C^{n) 
Jn Jn d 



and some positive constant C . Then there exists a positive constant C such that 

j (^|Vn|2 - dx>C J ^dx, Vug C^{n) . 

We finally mention the following related new Hardy inequality, which we think is of independent 
interest 

Theorem 3.10 For N > 3, let Q, C be a smooth bounded domain containing the origin, and define 
for e > 0, 

i^4^ if {x e n : d(x) > e} 

lid^y if {x^n-.dix) <e} . 
Then there exists eq = eo(17) such that for all < e < eQ and u £ Cq°{0,), we have 

Wul^dx > / VJxWdx . 



Proof: We will only sketch it. Let 0,i = {x £ : d{x) > e}. Then using the change of variable 

2-JV 

n := |x| 2 V, one can prove the following inequality 

,^ ,2 {N-2f 2\ r 2-N f 

- .1 19 ^ dx > -—- / --^x ■ V dS^. 

J 2 Jg^^ \x\^ 
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Similarly using the change of variable u := {x)X 2 (^d{x))v with X{t) = (1 — Int) ^ one can prove the 
following inequality 

1 \ 1 /■ 

V^l"" - , ,9, . u^ ]dx> - - / ^TT^c? • u dSx, 



n\ni V 4d2(x) / 4 Jqq^ d{x) 

for any < e < min{e~^,ei} where ei > is such that d^^X{d) + 2Ad(ln(i) > for d < ei. The result 

v > since 



then follows showing that for < e < eo = min{e ^, ei, 2N-3 } have ^^|jrp^x — ^i^Vd 

:= — Vd on dQi; here i? denotes a positive constant such that B{0,R) C 0,, which exists due to the 
assumption on // 



3.2 Complete sharp description of the heat kernel for small values of time 

In this section we prove the two-sided sharp estimate on the heat kernel k{t, x, y) stated for small time in 
Theorem 11.11 

Proof of Theorem li.il Since for any x £ 0, and for some positive constants ci, C2 we have the following 
estimate ci\x\^ d^ (x) < (pi{x) < C2\x\^ d'^ (x) for a = 2 and A = 2 — N, we can apply the result of Theorem 
12. 151 to the operator K = — :^^div{ipi{x)V). Hence due to 1)3. 11(1 the result follows immediately. 

II 

Let us finally make some remarks concerning Schrodinger operators having potential V{x) = c|x| ^. 
Arguing as in Lemma 7 in |DDj one can easily prove that the first Dirichlet eigenfunction for the 
Schrodinger operator —A — -j^, < c < -^^ — j^-, behaves like |x| 2d(x) on all where A := 2 — A*" + 
■sj {N — 2)2 — 4c. Then we have 

Theorem 3.11 For N > 3, let Q. C be a smooth bounded domain containing the origin. Then there 
exist positive constants Ci, C2, with Ci < C2, and T > depending on 0, such that 

Cimin|(|x| + + Vi)^'^ , ^^^^M^||a;y|lt-f e-<^2^^^ < 

, / N ^ . f/i I ^^ 1^1 /I I l-^l d(x)d(v)^, ,A JV ^, \x-y\^ 

< kc{t,x,y) < Cammjdxl + + Vt)^,^^^^||a;y|2i~^e~^i^^ , 

for all x,y £ n and < t < T ; here kc{t, x, y) denotes the heat kernel associated to the operator —A — 
in under Dirichlet boundary conditions for < c < -^^ — 4"; \ := 2 — N + {N — 2)^ — Ac. 

Theorem 3.12 For N > 3, let Q C be a smooth bounded domain containing the origin. Then there 
exist two positive constants Ci,C2, with Ci < C2, such that 

Ci d{x) d{y) |xy|te-^i* < k^{t,x,y) < C2 d{x) d{y) \xy\U~^^* , 

for all x,y £ ^} and t > large enough; here kc{t,x,y) denotes the heat kernel associated to the operator 

'-2 
4 



/- 2)2 

—A — in $7 under Dirichlet boundary conditions for < c < -^^ — 4"; '^i its (positive) elliptic first 
eigenvalue and X := 2 — N + ^ [N — 2)^ — 4c. 
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4 Critical boundary singularity 



In this section we prove Theorems 11.31 and 11.41 as weh as a new Hardy- Moser inequaUty (Theorem I4.3|l . 
The structure of this section is as follows. 

In Subsection 4.1 we first prove the improved Hardy-Moser inequality. Then in Subsection 4.2 we get 
the global in time pointwise upper bound for the heal kernel of the Schrodinger operator — A — (l/4)d~^(x), 
which is sharp when x and y are close to the boundary (see Theorem 14.4(1 . Then arguing as in |Dlj . we 
deduce the sharp heat kernel lower bound for time large enough, thus completing the proof of Theorem 

The proof of Theorem ll.3l is finally completed in Subsection 4.3, using the parabolic Harnack inequality 
up to the boundary stated in Theorem 11.51 



4.1 The improved Hardy-Moser inequality 

Here we will prove a new improved Hardy-Moser inequality which we think it is of independent interest. 
The proof is based on an auxiliary Hardy-Sobolev inequality, that we will show here, as well as on the 
following improved Hardy inequality stated in Theorem A in |BFTlj . 

Proposition 4.1 For N > 2, let Q C be a smooth bounded and convex domain. Then there exists 
Dq positive such that for all D > Dq 

2 f d{x) 



^ ,2 1 2\ , I f ^ \~J 2 , 

~ , ,9/ \ U dx> - — -r— — u dx , 
for any u G C^(0); here X{t) := t £ (0, 1]. 

Let us now state the auxiliary Hardy-Sobolev inequality we will use in the sequel. 

Lemma 4.2 Let a > 0, N > 2 and C be a smooth bounded domain. Then there exist i5o > and 
C = C{a, 6o) > such that 

N-l 

r , / r aN N \ N 

d"(x)|Vi;| dx+ / d''-^{x)\v\ dx > C [ / d—{x)\v\— dx) 

for any v G C^{^1) and any < 6 < 6o; here Jl^ := {x G $7 : dist(2;,3il) < 6}. 

Proof: We will follow closely the argument of |FMT2j . Our starting point is the following Gagliardo- 
Nirenberg inequality (see p. 189 in [Mj) 

^Tvii/ii^^^^^ < iiv/ii^i(f,) , V / G c^m , 

where Sn is a positive constant depending only on N. 

For any v G C^{Q) let us apply the above inequality to the function / := d'^v. Hence we obtain 

ra— 1 / 

a 

n 



SN\\d''v\\^^^^^< I d''{x)\Vv\dx + a I d''~'{x)\v\dx . 



We next estimate the last term above. Let ips G 0^(0,25), < 995 < 1, be a cut off function which is 
identically one in Qs and identically zero in R^ \ il25- Clearly v = ipsv + (1 — ips)v. Then we have 

Q / d''-^{x)\v\dx <a [ d'^-'^{x)\^sv\dx + a [ d''-^{x){l - ips)\v\dx < 
Jn Jn Jn 
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< a / ^ {x)\ipsv\dx + a / ^{x)\v\dx . 

Jn Jn\ns 

Concerning the first term on the right hand side we have 

a / d°'~'^{x)\(psv\dx = / Vd"" -Vdlifsvldx = - / d°'{x)\7d ■ \7\fsv\dx - / d°'{x)Ad\tfsv\dx < 
Jn Jn Jn Jn 

< [ d''{x)\V{ipsv)\dx + co5 [ d''-^{x)\^sv\dx , 
Jn Jn 

here we used the smoothness assumption on $7 which imphes that \dAd\ < cqS in Qs for 5 small, say 
< S < 5o, and for some positive constant cq independent of 5 {5o, cq depending on Q). Thus we have for 
any < S < 6o 

a [ d'^-\x)\ipsv\dx < ^— / d''{x)\V{ipsv)\dx <C [ d'^{x)\Vv\ipsdx+ 

Jn " - cqOo Jq 

+^ / d'^{x)\v\dx <C [ d''{x)\Vv\dx + C [ d'^-\x)\v\dx , 
Jn2s\ns Jn Jo2s\ns 

from which the result follows. 

// 

We next state the new improved Hardy-Moser inequality. 

Theorem 4.3 (Improved Hardy-Moser inequality) For N > 2, let Q, C be a smooth bounded 
and convex domain. Then there exists a positive constant C such that 

2_ 

jy (^|Vnp - -^u^^ dx^ u'^dx^ - ^ j ^^^^^^^^^ ' ^ ^ ^o'i^) ■ 
Proof: Changing variables hy v := ud~2^ we get 

2(^^+2) f JV+2 2(JV+2) f aN_ 2a\-^ 7 

u N dx = d N V N dx = d^-^{v )^~^dx , 
Jn Jn 

with a := (^+^K^~^) ^ Applying Lemma 14.21 to the function we have 



/ d^{v'^°')^dx <C I I d"|Vt;2"|(ix+ / d""^ 
Jn \Jn Jn\ns 



N 
N-1 



v^'^dx I < 

n\ns 



<C[2a / d"\Vv\\v\^''-^dx+ / d^'^v^^dx < 
\ Jn Jn\ns J 



N 
N-1 



N 



[ d{a 
Jn 



1 . . 1 / 2 \ 5 / \ ^ 



'2 I I „ -, „/„ ,N \ 2 



n / \Jn\ns d j \Jn\ns 



<Cl( i d{x)\Vv? dxV [ I d^^~^v^^^^-^UxY + { i ^dxl (/ d2"-i^;2(2"-i)dx 1 \ < 
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here we used Proposition 14.11 (observe that jX'^{j^) > jX'^{jj) if x G \ ils) and standard estimates. 
Returning to the original variable u, we obtain 

JV JV 



that is, 

2(JV-1) 



If = 2 we have that a = 1, thus the above inequality becomes 
which is the sought for estimate. For > 3, we use Holder inequality to obtain 

2(JV-1) _2_ JV-2 

2(iV+2) X'V / f \ ^ / f 2(N+2) \ ^ ( f / , ,9 1 o , 



from which 



2(iV+2) ^ , ^ I I 2 ' 1 ^ 



and this completes the proof of Theorem 14.31 



2 



u^^dx <C\ I u'dx ) { I [ \Vu\' - j^u' ) dx 



II 



4.2 Boundary upper bounds and complete sharp description of the heat kernel for 
large values of time 

Here we will first prove the following: 

Theorem 4.4 For N > 2, let Q C be a smooth bounded and convex domain. Then there exists a 
positive constant C such that 

t2 

To this end we need the following estimate of |i:''MT2j : 

Proposition 4.5 For N > 2, let Q, C be a smooth bounded and convex domain. Then there exists a 
positive constant C such that 

2 

(iV^^ - 4^^') dx>C (|^d§(^-2)-^(x)|^r dxy , (4.1) 

for any u E CQ°(r2) and any 2 < q < if N > 3 or any 2 < q < oo if N = 2. 
Using 1)4.11) the following logarithmic Sobolev inequality can be easily obtained 
I V \ f I. .0 I o\ ( iV+1 



v^\og\^^\jddx<e j \^vYd- -l\dv^\dx^\Kx — loge ) , (4.2) 
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for all V E C^(ri), V > 0, and any e > 0; here Ki is a positive constant independent of e and ||f||2 := 
To obtain we apply (|I?T|) to u := ^^"2 to get for any v £ C^(f^) 



2 



Taking g := '^'^^^l'* we have 



JV-l 



y [\Vv\'^d - ^Ad dx>c(^j v^'i^-') ddx^ ""^^ . (4.3) 
Then arguing in a quite similar way as in the proof of 1)3 .7(1 in Subsection 3.1 we obtain (|4.2() with 
Proof of Theore'm \4-4\ ' Let Hq{Q, d dx) be the closure of CQ°(ri) with respect to the norm 

Let H := U^^HU, U : L'^{Q,d dx) L'^{Q) being the unitary operator Uv := d^v, thus H := 
— ^div{d'V) — To the operator H defined in the domain D{H) = {f G HQ{0,,d dx) : Hv G 

L'^{Q,d dx)} it is naturally associated the bilinear symmetric form defined as follows 'H[vi,V2] ■=< 
Hvi,V2 >L2{n,d dx)=< vi,V2 >H^(n.d dx) which is a Dirichlet form. Whence Lemma 1.3.4 and Theorems 
1.3.2 and 1.3.3 in |D4j implies that e~^^ , which is an analytic contraction semigroup in L'^{Q., d dx), is also 
positivity preserving and a contraction semigroup in LP{Vl,d dx) for any 1 < p < cxd. As a consequence 
for any t > and any p > 2 

e-^*[L^{n,d dx) n L^{n)]+ c [H^{n,d dx)nLP{n,d dx)nL'^{n)]+ ; 

thus by density argument the logarithmic Sobolev inequality, which can be deduced as usual from 
the logarithmic Sobolev inequality ()4.2p (see Subsection 3.1 where a similar argument is used) more 
generally applies to any function in U(>o e~^*[L^(0,d dx) nL°°{Q,)]'^. This means that Theorem 2.2.7 in 
|D4j can be applied, as in Corollary 2.2.8 in |D4j . to the operator H; whence obtaining that 



and by duality that 
that is 

Here we use the following notation: 



0+ N+l 

e-^'^-^oo < Ct-— 



|e-^*||i^2<Ct-^ 

0+,, N+l 



-Ht\ 



\e \\q-^p-= sup 



0<11/|L<1 MM-^^Mg 



e-^7(^)ll; 
\\f{x)\\, 



where := (/f^l/l'^ ddx)''. This implies, by Dunford-Pettis theorem, that the semigroup e ^* is 

indeed a semigroup of integral operators; that is a heat kernel h{t, x, y) associated to the semigroup e~^'' 
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— 1 iV 

is well defined and satisfies tfie following pointwise upper bound h{t,x,y) < C^t^^, for any x,y Q 

and any t > 0. Theorem 14.41 then follows, due to the fact that, as a consequence of the unitary operator 
U, the heat kernels h{t,x,y) and h{t,x,y), corresponding respectively to H and H, satisfy the following 
equivalence h(t,x,y) = {x)d2 (^y) h(t,x,y) . 

Remark 4.6 Applying Davies's method of exponential perturbation to the operator H (see Section 2 in 
for details) the upper bound in Theorem \4-4\ can be improved by adding a factor cse ^ 

// 

Let us now give the sketch of the proof of Theorem 11.41 

Proof of Theorem \1.4[ Let us first improve by an exponential decreasing in time factor the global 
in time upper bound stated in Theorem 14.41 To this end let us define, as in Section 2 of |D2j . the 
operator H := U~^{H — Xi)U, U : L'^{Q,ip1dx) — > L'^{Q) being the unitary operator Uw := (piw, 
thus H := — ^dzf((/9?V). Here > denotes the first eigenfunction and Ai > the first eigenvalue 

corresponding to the Dirichlet problem —Aipi — ^^^^^-^ y^i = ^ifi in Q, ipi = on dQ, normalized in such 

a way that Vi{x) dx = 1. Since it is known that there exist two positive constants ci, C2 such that 

ci d^x) < (^i(x) < C2 V X G (4.4) 

(as a consequence of Lemma 7 in |DD) ). logarithmic Sobolev inequalities analogous to 1)4. 2|) also hold 
true if we replace H hy H. Thus as a consequence of Gross theorem, exactly as in the proof of Theorem 
14.41 the corresponding heat kernel h{t, x, y) satisfies the same pointwise upper bound as h{t, x, y) that is 
h{t, X, y) < ^-t~~ for any x, y G O and any t > 0. Prom the definition of U , it follows 

h{t,x,y) = ipi{x)ipi{y)h{t,x,y)e-^^* , (4.5) 

thus we get that h(t,x,y) < C "^^ (^)^^ iv) ^- — ^-Xit £qj. ^^^^y x,y £ Q and any t > 0. Finally arguing as in 

Theorem 6 of |Dlj , an analogous lower estimate can be easily deduced (we refer to the proof of Theorem 
11. 21 where a similar argument is used). 

// 

4.3 Complete sharp description of the heat kernel for small values of time 

In this section we prove the two-sided sharp estimate on the heat kernel h{t, x, y) stated for small time in 
Theorem 11.^^1 

Before doing so let observe that Theorem 11.51 entails also the following parabolic Harnack inequality 
for the Schrodinger operator having critical singularity at the boundary 

Theorem 4.7 For N > 2, let Q, C be a smooth bounded and convex domain. Then there exist 
positive constants Ch o,nd R = R{Vt) such that for x€Q,0<r<R and for any positive solution u{y, t) 
'^^ W ~ ~^ 4dF(y)^ {^{^7 r) n $7} X (0, r^) we have the estimate 

'^^(y,t)6{e(x,r)nc}x(4,^) u{y,t)d-'Hy)<CH ess ^n/(^ .^^^^(^ .)nn}x(|r-2,r-2) u{y,t)d-h (y) . 
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Proof: As a first step let us observe that if u satisfies Ut = —Hu then v{y,t) := e^^^ipi{y)~^u{y,t) 
satisfies vt = —Hv. Whence as a consequence of 1)4. 4() . due to Remark 12.161 v satisfies Theorem 11.51 for 
a = 1. Prom this the claim can be easily deduced. 

// 

The proof of Corollary 11.71 is similar to the above proof of Theorem 14.71 thus we omit the details. 
We are now ready to prove Theorem 11.81 

Proof of Theorem Since for any x £ 0, and for some positive constants ci, C2 we have the following 
estimate cid^^x) < V'i(x) < C2d2(x) for q = 1, we can apply the result of Theorem 12. 151 to the operator 
H = — :^j-^div{ipi{x)'V). Hence due to ()4.5() the result follows immediately. 

// 

The proof of Corollary 11.81 is similar to the above proofs of Theorems 11.31 and 11.41 thus we omit the 
details. 

Let us finally make some remarks concerning Schrodinger operators having potential V{x) = cd^'^(x). 
Arguing as in Lemma 7 in |DDj one can easily prove that the first Dirichlet eigenfunction for the 
Schrodinger operator —A — ^|3j^, < c < j, behaves like d^ on all for a := 1 + Vl — 4c. Then 
we have: 

Theorem 4.8 For N > 2, let ft C be a smooth bounded and convex domain. Then there exist positive 
constants Ci,C2, with Ci < C2, and T > depending on Q such that 

C, min 1 1, dHx)^ | < < | ^ dHx)^{y) ]^ ^-f ^-c.^ 

for all x,y £ 0, and < t < T; where hc{t, x, y) denotes the heat kernel associated to the operator 
—A — in under Dirichlet boundary conditions, for any < c < | and a := 1 + ^/l — 4c. 

Theorem 4.9 For N > 2, let C tl^ be a smooth bounded and convex domain. Then there exist two 
positive constants Ci,C2, with Ci < C2, such that 

Ci (it(x) dt(y) e'^i* < hc{t,x,y) < C2 d'^ix) d"^ (y) e'^i* 

for all x,y £ ^} and t > large enough; where hc{t, x, y) denotes the heat kernel associated to the operator 
—A — ^^^^ in Q under Dirichlet boundary conditions, Ai its (positive) elliptic first eigenvalue, for any 

< c < i and a := 1 + y/l - 4c. 

Remark 4.10 Let us at this point remark that Theorems \l.cl\ and \4.f^ as well as Theorem \4. 7| concerning 
respectively sharp asymptotic for small time and the parabolic Harnack inequality up to the boundary for 
the Schrodinger operator having potential V{x) = cd~'^{x), hold true also without any convexity assumption 
on the domain Q under consideration. 



4.4 On Davies conjecture 

In this subsection we consider Davies conjecture. For this we suppose that E denotes the self-adjoint 
operator associated with the closure of the positive quadratic form 

i,j=i J / 
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initially defined on C^{0,); where {aij{x))nxn is a measurable symmetric uniformly elliptic matrix such 
that 

N 

«J=1 



and y is a potential on 0, such that 



where 



V{x)=Vi{x) + V2{x) , (4.1) 



l^i(x)l<^^, V2{x)eLP{n), p>^ . (4.2) 
We also suppose that 

Ai := inf > 0, (4.3) 

o^v3GC^(n) ip^dx 

and that to Ai there corresponds a positive eigenfunction ipi satisfying for all x € ^2 the following estimate, 

ci d^{x) < ipi{x) < C2 d^{x) , for some a > 1 (4-4) 

and for ci, C2 two positive constants. 

Thus E is defined on the closure of Cq°(0) with respect to the norm defined by the quadratic form Q. 
Then as before it can be shown that is a well defined nonnegative self-adjoint operator on L'^{Q) such 
that for every t > 0, e~^* has a integral kernel, that is e~^^uo{x) := f^e{t, x, y)uo{y)dy and if iV > 3 a 
Green function Gj^{x,y) = J^e{t,x,y)dt denoting the kernel of E~^. 

Theorem 4.11 For N > 3, let ^ C he a smooth hounded domain. Suppose that i4-^) , i4-3\) 

and i4.4}) are satisfied. Then there exist two positive constants Ci,C2, with Ci < C2, such that for any 
x,y £ Q, 

^ . f 1 d^{x)d'^{y) 1 ^ , , ^ . f 1 dt(x)(it(y) 1 

Ci mm < ■; , „ , ■; n^T- > < G^(X,V) < Go mm < ■; , „ , ■; n^T- ;r > . 

Davies conjecture is stated under slightly stronger assumptions on V than ()4.2() and on ipi (only when 
a = 2). 

Proof: We note that we have Ei := — l^div{ipjV) = U-'^{E - \i)U, U : L'^{Q, ipfdx) L'^{Q) being 
the unitary operator Uw := (fiw; hence we have the following relationship between heat kernels 

e{t,x,y) = ipi{x)ipi{y)ei{t,x,y)e~^^^ . (4.5) 

Due to ()4.4I) we can apply Theorem l2.15l to the operator Ei. Hence due to 1)4. 5|) for two positive constants 
Ci ^ C2, we have for small time 



d2 {x)d2{y) \ I d2{x)d2{y) \ N 



2 



Cimin |1, — ' '^^ I t-^e-^^—r- < e{t,x,y) < C2min|l, — ' t'^e-^^—^. (4.6) 

On the other hand for large time 

Ci (it(x) dt(y) e~^i* < e{t,x,y) < C2 dt(x) (it(y) e'^i*, (4.7) 

for all x,y G 0,. To obtain this estimate we need to prove a global Sobolev inequality on fi, which can be 
easily deduced from its local version 1)2.12(1 as well as (|2.18|) with A = there, by means of a partition of 
unity as in [K]. Then the result follows integrating e{t,x,y) in the time variable. 
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